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An analysis method for transmission measurements of superconducting 
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Superconducting resonators provide a convenient way to measure loss tangents of various dielectrics at low 
temperature. For the purpose of examining the microscopic loss mechanisms in dielectrics, precise measure- 
ments of the internal quality factor at different values of energy stored in the resonators are required. Here, 
we present a consistent method to analyze a LC superconducting resonator coupled to a transmission line. 
We first derive an approximate expression for the transmission S-parameter 6*21(0;), with uj the excitation 
frequency, based on a complete circuit model. In the weak coupling limit, we show that the internal quality 
factor is reliably determined by fitting the approximate form of 5'2i(w). Since the voltage V of the capacitor 
of the LC circuit is required to determine the energy stored in the resonator, we next calculate the relation 
between V and the forward propagating wave voltage V-^, with the latter being the parameter controlled in 
experiments. Due to the dependence of the quality factor on voltage, V is not simply proportional to V-^. 
We find a self-consistent way to determine the relation between V and V^^, which employs only the fitting 
parameters for 6*21(0;) and a linear scaling factor. We then examine the resonator transmission in the cases 
of port reflection and impedance mismatch. We find that resonator transmission asymmetry is primarily due 
to the reflection from discontinuity in transmission lines. We show that our analysis method to extract the 
internal quality factor is robust in the non-ideal cases above. Finally, we show that the analysis method on 
LC resonator can be generalize to arbitrary weakly coupled lumped and distributed resonators. The gen- 
eralization uses a systematic approximation on the response function based on the pole and zero which are 
the closest to the resonance frequency. This Closest Pole and Zero Method (CPZM) is a valuable tool for 
analyzing physical measurements of high-Q resonators. 
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I. INTRODUCTION 

Measuring the response of a resonator is a generic 
method for measurements of properties of solids at low 
temperature, such as surface resistance^ of superconduc- 
tors and complex permittivity^ of dielectrics. It is also 
very useful for sensitive detection of photons'^, states of 
superconducting quantum bits'*, and the zero-point mo- 
tion of mechanical oscillators^. Recently, the need for in- 
creasing coherence times of superconducting qubits mo- 
tivated new interest in searching for low-loss dielectric 
materials as well as understanding their loss mechanism. 
Internal quality factors of several types of superconduct- 
ing thin-film resonators^*^ have been measured in detail. 
The most common way to measure the internal quality 
factor is to measure the S-parameters of resonant circuits. 
For this purpose, a number of methods were developed to 
extract the quality factor of a single resonator from reflec- 
tion or transmission measurements*^. In order to mea- 
sure multiple resonators simultaneously, methods*'*^'** 
based on frequency multiplexing have been developed. 
However, all the previous methods derive the scatter- 
ing parameters (S-parameters)*^ by using equivalent cou- 
pled circuits in the vicinity of the resonance. Because 
the equivalent circuit is not uniquely defined, different 
approximations involving such circuits lead to different 
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reduced S-parameters. This results in different expres- 
sions of the internal quality factor. In this paper, we 
present a consistent method to approximate the exact S- 
parameters by analyzing the algebraic properties of the 
S-parameter function. We start from a circuit model in 
which a LC resonator couples to a transmission line both 
inductively and capacitively. We develop the Closest Pole 
and Zero Method (CPZM) to reduce the transmission S- 
parameter 6*21(0;), with o; the excitation frequency, of 
the circuit to a ratio of complex linear functions. The 
approximated 6*21(0;) only depends on four parameters: 
resonance frequency ojq , internal quality factor Qi , exter- 
nal quality factor Qe and a parameter Qq which is related 
to the resonance asymmetry. The internal quality factor 
Qi can be extracted from fitting the data to the simplified 
6*21(0;) function as opposed to the full circuit-dependent 
form. 

We also discuss another important aspect of dielectric 
loss characterization. In order to understand sources and 
mechanisms of the loss, the dependence of dielectric loss 
on energy stored in the resonator is required. The energy 
dependence of loss can be characterized^'*°'^° by measur- 
ing Qi as a function of the voltage V across the capacitor 
in the LC resonator. However, this is a nontrivial task 
since one usually does not have precise knowledge of all 
the circuit parameters needed to calculate V. We show 
that the ratio V/V-^, where V-^ is the forward propagat- 
ing wave amplitude at the excitation port, can be cal- 
culated using parameters determined from the fit of the 
6*21(0;) function and an additional scaling factor. The 



scaling factor depends on circuit parameters which can- 
not be directly determined from the fit. However, it can 
be estimated based on circuit parameter simulations. An 
important point to note is that the knowledge of the volt- 
age V^ up to a fixed factor is still valuable for comparing 
experimental results with theoretical predictions of the 
loss mechanism. 

Although the reduced S21 fits the experimental data 
extremely well, in some measurements the resonance 
curve shows a pronounced asymmetry, which is incon- 
sistent with realistic values of circuit parameters. We 
analyze a more complete model which takes the reflection 
of wire bonds and impedance mismatch into account. We 
find that for realistic circuits the asymmetric resonance 
primarily arises from wave reflection by the wire bonds 
used to connect the device and the measurement appa- 
ratus. We show that despite these non-idealities, our 
procedure for extracting Qi based on the simplified 6*21 
form is still valid. 

We note that, in previous work, simplified equivalent 
circuits were used to obtain a simple expression for the 
S-parameters in the linear fractional form^^. However, 
in this paper, we develop a systematic approach to re- 
duce the response function of an arbitrary resonator to 
a linear factional form by using specifically its algebraic 
properties. 

This paper is organized as follows. In Section H, we 
approximate the transmission S-parameter 6*21 of the res- 
onance circuit we consider to a linear fractional form by 
applying the CPZM. In Section HI, we investigate the 
complex function v{uj) — V/V(^^, where V is the capaci- 
tor voltage and T^^ is the forward propagating wave volt- 
age. We approximate v{uj) to a linear fractional form 
using the same method and we find that v depends on 
Qi in general. We derive a formula to transform the 
points {V(^,Qi) to {V, Qi). In Section IV, we employ the 
CPZM to analyze the transmission parameter 6*21 of a 
circuit where we consider the inductance of the bonding 
wires and the effect of variation of impedance of different 
transmission lines. In Section V, we discuss the gen- 
eral application of CPZM to the circuit response of both 
lumped and distributed high-Q resonators. 



II. DERIVATION OF RESONATOR TRANSMISSION 



transmission lines have characteristic impedance Zq. 



B. Exact 521 from circuit model 

Fig. 2 shows the equivalent circuit of the network 
shown in Fig. 1(b). Since the impedance of the source 
and load match the characteristic impedance of the trans- 



mission line, we have Va = 21^„ and VI, 



Kut where 



Vg is the voltage of the source and T4ut is the voltage 
on the output load. We use Kirchhoff 's laws to calculate 
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FIG. 1. (a) Optical microscope image of an LC resonator 
coupled to a transmission line, designed for measuring the 
loss tangent of the dielectric in the parallel-plate capacitor. 
The device is fabricated using two aluminum layers, patterned 
with electron-beam lithography. The dielectric is a 5 nm thick 
layer of AlOa;. (b) Lumped element circuit model for the 
device in (a). 



A. Experimental design and circuit model 

We consider a LC resonator coupled to a transmission 
line. In our experiments, the resonator consists of a 
meander inductor L and a parallel capacitor Cq. The 
transmission line is a coplanar waveguide (CPW) (see 
Fig. 1(a) for the sample design). We use a model 
circuit similar to that proposed in Ref. 18, as shown in 
Fig. 1(b). The LC resonator couples to the transmission 
line through a mutual inductance M and a capacitor 
Cc- Dielectric loss is modeled by a resistor R shunting 
the lumped capacitor Cq. Both the input and output 
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FIG. 2. Equivalent circuit of the circuit showed in Fig. 1(b). 



^21 - V+^JV+ = 2Kut/^5, which gives: 
^21 ^ f 2 + il (juoC,Z[^ + ^ 



(1) 



where 



and 



V 



_2f -jc^CZ^t 




Z[^ = Zo+:]ojLi-]ujAP/L, 
7' = -^0 

°"* 1+JC.CeZo' 



and Z(7 = 



1 



jujCa + l/R' 
C. Reduction of >S'2i to a linear fractional form 



(3) 
(4) 

(5) 



TABLE I. Parameters of the circuit shown in Fig. 2 



Li 



M 



a 



Co Zq R 



0.71 pH 11.9 pH 288.7 pH 5.0 fF 2.5 pF 50 » 50 kP, 

TABLE IL Reduced parameters (see main text for explana- 
tions). 



a 



P 



7 



0.041 0.0025 0.0020 0.21 0.00021 



The S-parameter ^21 in Eq. (1) is in a polynomial frac- 
tional form of uj. Our goal is to derive, in a systematic 
manner, a simpler approximate form of the 6*21 . Since 6*21 
is dimensionless, it is convenient to express it in terms of 
a set of dimensionless parameters. We choose the follow- 
ing dimensionless parameters: 
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(6) 



Here Wq = l/\/LCo is the resonance frequency of the 
uncoupled resonator. Coupling to the transmission line 
will result in a shift in the resonance frequency from Wg. 
When expressed in terms of the new parameters, the ex- 
pression for 521 becomes: 



^21 



/(^) 



flo + aiqx + a2{qxy 



g{x) bo + biqx + b2{qx)^ + b3{qx)^ + b4{qxy 

(7) 
where f{x) and g{x) arc polynomial functions of qx. The 
expressions for the coefficients in Eq. (7) are given in 
Appendix A. 

The dimensionless parameters in Eq. (6) arc chosen 
such that their values are much smaller (i.e. a, /3, 7, 
and q) or comparable (i.e. x and ^) to unity, which 
allows for developing systematic order expansions. To 
illustrate this, in Table I, we list realistic values of 
circuit parameters for a resonator realized in our ex- 
periments with an estimated bare resonance frequency 
w'q = 2tt X 5.9 GHz. Corresponding values of the reduced 



parameters are shown in Table IL We note that while all 
the reduced parameters arc small, q may be smaller than 
one by many orders of magnitude (the quality factor of 
a superconducting resonator often reaches 10^ especially 
when driven at high power). For resonators weakly cou- 
pled to the transmission line, we expect that the reso- 
nance frequency and the width of the resonance are not 
perturbed very strongly. For this reason, the relevant 

range of the variable x = 2Q " ,"" corresponds to \x\ < 1. 

Since qx is a number much smaller than the rest of 

the parameters, we can approximate f{x) and g{x) by 

keeping only the terms up to the first order of ga;, which 



gives: 



f{x) _ ao + aiqx 
g{x) 60 + biqx ' 



(8) 



Fig. 3 illustrates the accuracy of the approximation above 
by comparing both the amplitude and phase of 521 and 
S21 using realistic circuit parameters given in Table I. 
Noting that both f{x) and g{x) contain terms up to the 
first order in x, we rearrange the terms and write S21 in 
the following linear fractional form: 
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where a= ^^ 

q ai 



Aole^'i' = ai/bi. We note that 



= ill' l^oie- = ui/oi. 
the parameter |Ao|e-''^ has a value which is very close to 
1 . For this reason it would be difficult to distinguish this 
factor from attenuation and gain factors in the measure- 
ment setup. Therefore, we simply drop this factor from 
now on. The expressions for a and b are ratios of poly- 
nomials, depending on all the reduced parameters, given 
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FIG. 3. Amplitude (a) and phase (b) of the exact S21 (red 
soUd line) and the reduced linear fractional form 521 (blue 
dashed line). The expressions are defined in Eqs. (7) and (8), 
respectively. 



in Appendix A. We rewrite next the parameters a and b 
in a way which makes the dependence on q explicit. We 
have: 



and 



, B 



D F 

b = C+-+jiE + -), 

q q 



(10) 



(11) 



where A, B, C, D, E and F depend on the reduced 
parameters other than q. 

We simplify the coefficients A, C, and E by using a 
Taylor expansion in terms of a and 7, which characterizes 
the resonator coupling and are therefore very small. We 
obtain: 



A =1 + 0(7), 

C=l + 0(7)+0(a2), 

E=0{-i) + 0{a^). 



(12) 
(13) 
(14) 



This expansion is carried out using Mathematica'^^ . We 
only retain the zero order terms for A, C, E. A first 
order term for E is not needed, due to the large F/q 
term. With these approximations, the function 6*21 takes 
the form: 



'S'21.rcd — 



J q 



JX 



1 + -+J- 
q -J q 



JX 



(15) 



The resonance curve described by Eq. (15) is asym- 
metric in general. We identify the resonance frequency 
as the value when |S'2i.rcd| is minimized. The resonance 
frequency wq depends on the internal loss of the capaci- 
tor (details shown in Appendix A). However, for high-Q 
resonators, this dependence is weak for (7 <C 0. For sim- 
plicity, wc take the resonance frequency as the resonance 
frequency without internal loss (q ~ 0) as the following: 

uj'oil-B/2), (16) 



where 



Wo 
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7(1 -a) 



1 + 7 — a7 
The internal quality factor is defined as 
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i{l-B/2)RCo 



(17) 

(18) 
(19) 



We rewrite Eq. (15), showing specifically the resonance 
frequency cjq and the internal quality factor Qi to 



1 + 2jQ, 



'21, red — 



1 



§+j-^ + 2jg. 



(20) 



with Qa and Qe defined as the following expressions: 

(^a - ^ _ ^ , 

1 - -^/2 



D 



(21) 
(22) 



The reduced S21 m Eq. (20) is in a form which is similar 
to a Lorentzian. Qa is a parameter that characterizes 
the asymmetry of ^21 with respect to the resonance fre- 
quency. In the limit where Qa -> 00, ujo{Q~ + Qe^) 
is the width of the transmission single in frequency do- 
main. Therefore, Qe is the external quality factor of the 
resonator. 

Eq. (20) is the main result of this section. By fitting 
Expression (20) to 5*21 data, we can uniquely determine 
the fitting parameters wq, Qi, Qe, and Qa- It is also 
worth to point out that Eq. (13) in Ref. 18 and Eq. (3) 
in Rcf. 8 also possess similar linear fractional forms, but 
here we first develop a consistent approach on deriving 
this form. In Section V, we discuss the more general ar- 
guments that justify this systematic approach. We show 
there that this method approximates the S-paramcter 
from a complex function with multiple poles and zeros 
to a reduced form governed by only one pole and one 
zero whose real parts are the closest to the resonance fre- 
quency. We call this method the Closest Pole and Zero 
Method (CPZM). 



III. VOLTAGE ON LUMPED ELEMENT CAPACITOR 

Characterizing the field dependence of dielectric loss 
is essential in understanding the loss mechanism of di- 
electrics. For amorphous dielectric materials, tunneling 
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FIG. 4. Amplitude (a) and phase (b) of the exact V jV-^ 
(red solid line) and the reduced linear fractional form v (blue 
dotted line). 



two-level systems (TTLS)^*^'^"^ are believed to constitute 
the dominant source of loss in the quantum regime. The 
TTLS model has detailed predictions of the loss depen- 
dence on the AC electric field in the dielectric. In a LC 
superconducting resonator, the electric field in the dielec- 
tric is proportional to the voltage V on the capacitor. 

Using circuit theory, wc can determine ^(Qi,!^^), 
with Qi the resonator internal quality factor and V-^ the 
voltage of the forward propagating wave towards the res- 
onator input port. If Qi is a constant or the circuit is 
strongly coupled to the transmission line, V is propor- 
tional to V^- The experimentally measured (Q^,!/^^) 
dependence is transformed into a {Qi , V) dependence as- 
suming a constant ratio VjV-^^. In this section we develop 
a systematic method to express the dependence of V on 
V-^ in a way which takes into account the voltage de- 
pendence of the quality factor. We find that the (Q^, V^) 
dependence determined in this way is different than the 
case when V jV-^ is simply assumed constant. The dif- 
ference is more than a trivial scaling factor, and has an 
effect on the exponents appearing in the voltage depen- 
dence of the loss for realistic circuits. 

Starting from Eq. (2), we introduce v{lo) = VjV-^ and 
express it using the dimensionless parameters introduced 
in Eq. (6). In terms of the reduced parameters introduced 
earlier, u is a high order polynomial fractional form v ~ 
fv{x)/gv{x) (sec Eq. (A2) in Appendix A for the full 
expression of v). We apply our CPZM to approximate v 
to u, a linear fractional form of q: 



fvjx) _ cq + ciqx 
gv{x) do + diqx^ 



(23) 



with fv{x) and gv{x) containing terms up to first order 
in X. Fig. 4 illustrates the accuracy of the approximation 
above by comparing both the amplitude and phase of v 
and V using realistic circuit parameters given in Table I. 
We rewrite v as 



" = ^77- 



jx 



jx 



where C = ci/di, a' = ^^ and b' = ^^ 
from the rest of the reduced parameters gives 



(24) 



Separating q 



c 



K/q + jJ/q + jx 



H + N/q+j{0 + P/q)+jx' 



(25) 



Interestingly, some of the coefficients in the above equa- 
tion are exactly the same as those in the 6*21 formula 
given by Eq. (20). We have H = C, N = D, O = E, and 
P = F. Since it is the features of the near-resonance 
transmission that allow one to extract Qi, we calcu- 
late the value of v{uj) for a; = wq. At this frequency, 
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) = -|. Wc obtain: 


V 


A 


v+ 


,^ Vq ~~ 

,, ,, qi + qe+jqa 



K'^oM 



LJ—UJO 

where qi = l/Qi, qe = 1/Qe, qa = 1/Qa and 
.K + j{J-B) 



A = C- 



1 -B/2 



(26) 



(27) 



A is a scaling factor which does not depend on internal 
quality factor of the resonator. The full form of A is given 
by Eq. (A28) in Appendix A. 

To estimate this scaling factor, we write A as a Taylor 
series of small numbers 7, a, and /? to their first orders, 
which gives 



A = -a^ + J7 ■ 



(28) 



where S, has a quantity which can be comparable to 
unity. We approximate the expression above further to 
A « -Mio'JZo+jCc/Co « -Mloo/Zo+jCJCo. Given 
that loq can be obtained from fitting the S21 data and Zq 
is usually assumed to be 50 f2, we only need to estimate 
circuit parameters, M, Cc, and Co, for determining A. 

In experiments, one measures the S21 at various val- 
ues of V^J. For each value of V^^, a corresponding Qi is 
obtained from fitting Eq. (20) to the S21 data. Qe, Qa, 
and Wo are obtained from the fit as well. Then one can 
use Eq. (26) to transform data set (V^^, Q;) into {V, Qi). 
To illustrate the difference between the above two data 
sets, we plot in Fig. 5 the experimental (V^,^,(5i) data 
of the parallel plate LC resonator shown in Fig. 1(a) and 
the corresponding {V, Qi) data transformed with Eq. (26) 
with Qe = 1984, Qa = 4128, wo = 27r x 7.665 GHz, and 
A = 0.01146. The dielectric layer between the paral- 
leled plates is made of aluminum oxide, whose quantum 
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FIG. 5. Experimentally measured loss dependence on 
y;^(blue dots) and on y(red squares). The red line is a result 
of fitting Eq. (29) to {V, Qi) which yields A = -0.06f3 and 
Qifi = 397. The blue dashed line is a result of fitting Eq. (29) 
to (T/+, QO which yields A = -f .05 and 0,,o = 39f . 



regime loss is believed to be due to the TTLSs. Accord- 
ing to Rcf. 20, the TTLS loss at a constant temperature 
is predicted to be 



1 



l/Q^fl 



v/1 + iy/yc) 



2-A 



(29) 



with A = 0. We fit Eq. (29) to both (F, Q,) and {V+, Q,) 
data. We note that the fit to (F, Qi) yields A = -0.0613, 
which is a very good agreement with the TTLS model 
while the fit to {V-^^Qi) yields A = —1.05. As a result, 
the transformation from V-^^ to V with the loss of the 
resonator taken into account is crucial for determining 
the voltage dependence of dielectric loss. However, we 
find from Eq. (26) that V is approximately proportional 
to V-^ in the over-coupled limit, Qi » Qe- In this case, 
the fits to {V,Qi) and {V-^,Qi) will give approximately 
the same exponents A. 



IV. EXPLANATION OF ASYMMETRIC 
TRANSMISSION 

Although the function given in Eq. (20) fits the experi- 
mental data very well, the amount of transmission asym- 
metry characterized by Qa cannot always be explained 
by the circuit model in Fig. 1 (b) with a set of realistic pa- 
rameters. Fig. 6 shows two 15*211 curves, of LC resonators 
measured in our experiments, which are highly asymmet- 
ric. A similar, large, asymmetry was observed in experi- 
ments of Geerlings et al.^ and Megrant et al*. Khalil et 
al.^* attribute this effect to impedance mismatches. Our 
experiments indicate a very large impedance mismatch is 
required to explain the asymmetry transmission, which is 
unlikely to occur in our setup. Motivated by this obser- 
vation, we analyze a more complete circuit model, shown 



in Fig. 7. In this model we take into account the connec- 
tions between transmission lines on chip and outside the 
chip, implemented with bonding wires. We model these 
using inductances Lin and Lout- These connections act 
as reflection points for the propagating microwaves. The 
two parts of the on-chip transmission lines, between the 
bonding connection and the resonator, are assumed to 
have characteristic impedance 2^in/out and length /in/out 
respectively. 

For such a circuit, consisting of several two-port com- 
ponents, it is convenient to use the ABCD transmission 
matrix method^^. The ABCD matrix T of a general 2- 
port circuit is defined as 



(30) 





where VJ 



in/out 



and /; 



in/out 



are the voltages and currents 
at the input and output port respectively. We first use 
Kirchhoff's law to calculate the ABCD transmission ma- 
trix of the resonator part of the circuit (see Fig. 7 for the 
circuit schematics) . We write down the circuit equations 
which relate the voltages, Vi and V2, and the currents, 
/i and I2, to the capacitance voltage V as follows: 
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FIG. 6. Measured transmission curves (red dots) and the fits 
using Eq. (20) (black line) for two resonators on the same 
chip showing different types of asymmetries. The extracted 
fit parameters are indicated in the graph legends. 
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FIG. 7. Circuit model of a LC resonator coupled to a non-ideal transmission line. The on-chip CPW transmission line close 
to the input/output port has length ^in/out and characteristic impedance .^in/out- They are connected to the input/output port 
through wire bonds model by inductance I/in/out. Both the input and output ports are terminated with an impedance Zq. 



Based on these equations we can derive the ABCD matrix for the circuit 



where 



T., 




(34) 



Ta = 
Tb = 
Tc = 



-juj (L - CcLLiUj^ + CcAPuj'^) +R{-1 + {CqL + Cc{L + Li - 2M))uj^ + CqCc {-LLi + NP) cj*) 

-jLoj + R{-1 + {{Co + Cc)L - CcM)oj^) ' 

Lu {-jLiR + {LLi - RP) Lo + j{Co + Cc) {LLi - IVP) Rlu^) 
-jLuj + R{-1 + {{Co + Cc)L - CcMp^) ' 

CcLJ {Llu + jR (-1 + CoLio^)) 
-jLuj + R{-1 + {{Co + Cc)L - CcM)uj^) ' 

-jLuj + fl (~1 -f (Co + Cc)Luj^) 
-jLuj + R{^1 + {{Co + Cc)L - CcM)uj^) ' 
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For the other parts of the circuit, the ABCD matrixes of 
the wire bond inductors are 



(39) 



and the ABCD matrixes for the two pieces of the trans- 
mission lines between wire bonds and the resonator are 




TtLo 



cos Pal a jZa sin I3ala 

j sin 13 Ja/Za cos Pa L 



(40) 



where a can stand for "in" or "out", f3a = ^V LaCa, and 
La and Ca arc the characteristic inductance and capaci- 
tance of the transmission line. We define the time, for mi- 
crowave to propagate la in distance, as ta = laV LaCa- 
Given the width of the resonance Su <^ It: jta-, wc can 
approximate TtLq to an w-independcnt ABCD matrix 



-TL„ — 



COSWgtQ iZaSmUj'i^ta 

j sinUJQta/Za COSOJQta 



(41) 



I 

in the vicinity of a; = ( 

ABCD matrix of the complete circuit by multiplying all 

the ABCD matrices: 



In the end, we obtain the 



Lall 



TLi„TTLi„TTTLo 


,tTL„,t 


(42) 


TLi„TTLi„TTTLo 


utTio.f 


(43) 



The scattering matrix element 521 can be calculated from 
the ABCD matrix using the following equation: 



S" 



Tall, A + T^W^bIZo + TalKC-^0 + ^alLn 



(44) 



where Taii,a with a ^ A^ B^ C, or D is each of the matrix 
elements of Taii. 

In Fig. 8, we plot 5*21 in Eq. (44) using ijn = iout = 
370 pH, Zin — Zout = 50 ri, and parameters in Ta- 
ble I. These values of the inductance were determined 
based on numerical simulations of transmission through 
wire bonds with a 1 mm length. Fig. 8 shows that 5*21 
with different (?;„ = iin^o and ^out = ^out'^o '■^^^ have 
different amount of asymmetry with different signs. In 
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FIG. 8. The asymmetric resonance calculated using differ- 
ent values of transmission line electrical length tin/out • Here 
^in/out = i^o^in/out- Each curve has an incremental offset of 
—3 dB. The asymmetry degree parameter Qa is calculated 
for each curve. The symmetric 15*211 case corresponds to 



Qg and Q'^ here are different from the corresponding Qe 
and Qa in Eq. (20). Tfiey in general depend on more 
circuit parameters, i.e. L; 



in/out ■ 



^in/out a-nd tin/out- in 



addition to those circuit parameters in Table I. By fitting 
Eq. (46), the intrinsic quality factor Qi and the resonance 
frequency wq can be determined rehably. 

We emphasize that the same Qi and ujq appear in both 
Eq. (20) and Eq. (46) despite different external circuits 
coupled to the resonator. This shows that the impen- 
dence mismatch and wave reflection in the transmission 
line considered in this case does not compromise the de- 
termination of the internal quality factor Qi. 



V. GENERALIZATION TO OTHER RESONANCE 
CIRCUITS 

In this section we discuss the general application of 
CPZM to lumped and distributed resonators. We start 
with the consideration of lumped resonators. In this case, 
the response function (be it the transmission, reflection, 
or the ratio between voltages of any points) can be ex- 
pressed as a polynomial ratio: 



Rico) 



N{uj) 



(47) 



the multiplexed resonator measurement design, a num- 
ber of resonators are coupled to the same transmission 
line. Since their positions relative to the reflection points 
(wire bonds here) and their resonance frequencies are dif- 
ferent, the asymmetry feature of their transmission could 
also be different. This is what we observe in the experi- 
ment, which strengthens the hypothesis of our wire bond 
reflection model. The asymmetry is essentially due to 
reflections at the connection points which leads to weak 
standing waves with a length and frequency dependent 
pattern. 

Finally, we apply our approximation procedure shown 
in Section II to this more complicated model and simplify 
the transmission S21 in Eq. (44). We find that we reach 
a similar form to Eq. (15), which is given by 



''21, ro 



1+jf +ja^ 



' q -J q 



■JX 



(45) 



In the equation above, B is exactly the same as in 
Eq. (17). D' and F' are g-independent parameters. By 
using definition ujq = ui'q{1 — B /2) and Qi = looRCq which 
are the same as Eq. (16) and Eq. (18), and definitions 
g'g = (1 - B/2)/{F' - B) and Q'^ = il- B/2)/D' which 
are similar to Eq. (22) and (21), we obtain 



'-'21, red 



1 + 2jg, 






Q'e 






■2jQ.^^ 



(46) 



which has the same form as Eq. (20). The parameters 



where N{io) and D{io) are polynomials, of degree m and 
n respectively. We can factor both polynomials to obtain 



R{io) = RqX 



117=0' i^- Pi) 



(48) 



with Rq a complex number and Zi and pi zeros and poles 
of the response function R{uj). The entire resonance be- 
havior of the circuit is described by the positions of the 
zeros and poles and the complex scaling factor. Let us 
consider the response of the resonator close to one of 
its bare resonances, at frequency uj'q. If the resonator 
is weakly coupled to its measurement setups, we expect 
that the resonant behavior will be reflected by a zero 
and/or a pole of R{uj) which are close to lu'q. We denote 
the closest zero and pole by zq and po- The response 
function can be expressed in the following form; 



Riu;) 



N{uj){u^zo) 
D{uj){uj-po)'' 



(49) 



where N{uj) and D{uj) are polynomials with roots Zi^o 
and Pi^o respectively. For the near-resonance response 
of high-Q resonators, Nicj) and D(u;) arc approximately 
constants within range, ujq — Suj < u < ujq + Suj , where 
Sid = lo'q/Q. They can be well approximated with their 
values at cj = Wq. Therefore, the response function is 
reduced to 



i?(w) 



NH) 



zo) 



D{uj'o) (w-po)' 



(50) 
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FIG. 9. An asymmetric (a) and symmetric (b) response func- 
tion with its pole and zero plotted in the complex plane in (c) 
and (d) respectively. 



This approximate response function R{uj) only depends 
on the zero and pole zq and po in addition to a com- 



plex scaling factor 



D{u 



The response function R{uj) is 



asymmetric in general while a symmetric response func- 
tion requires ^[zq] = 5ft[po] (see Fig. 9). We note that 
any experimental measured frequency response, e.g. S- 
parameters, of a high-Q resonator can be fitted using 
Eq. (50) with three complex numbers, the zero, the pole 
and the complex scaling factor, as the fitting parameters. 
We also note that the approximation we applied in Sec- 
tion II, III and IV for reduced expressions for 521 and 
l^/Vjjj are equivalent to finding the approximate form of 
the closest pole po and zero zq in Eq. (50). 

For distributed resonators such as coplanar waveguide 
resonators, the response function is not simply a ratio of 
polynomials of w but also depends on a wave propaga- 
tion term z = e^^'^ , where r corresponds to the electrical 
length of the waveguide. The general response function 
can be written as 



R{u:) 



N{io,z) 
D{io,z)' 



(51) 



VI. CONCLUSION 

We developed the Closest Pole & Zero Method 
(CPZM) for analyzing the near-resonance response of su- 
perconducting resonators. We showed that in the high- 
Q limit, the response functions are well described by a 
linear fractional form. We first considered the transmis- 
sion S-parametcr (521) of a lumped resonator. Using the 
CPZM, we obtained a linear fractional 521 of this cir- 
cuit. By fitting the linear fractional ^21, one can extract 
the resonance frequency, internal quality factor, and ex- 
ternal quality factor out of the experimental measured 
transmission. We then applied the CPZM to analyze the 
relation between the capacitor voltage V and the exci- 
tation voltage V^^ of the same circuit. We found that 
the ratio between V and V^^ depends on all the fitting 
parameters of the transmission function in particular the 
internal quality factor. Since the internal quality factor 
depends on V itself due to the microscopic mechanism of 
the loss, V is not proportional to V^^ in general. We pre- 
sented a self-consistent method to calculate V from the 
fitting results to the transmission S-parameter for every 
given V^. 

Moreover, we studied the case in which imperfections 
of the transmission lines used for S parameter measure- 
ments are modeled by introducing impedance mismatch 
and inductance from wire bonds. We found that the 
asymmetric transmission near the resonance is mainly 
caused by wave reflection on the wire bond positions. 
Using the CPZM, we showed that our fitting routine for 
extracting the internal quality factor applies to this non- 
ideal case as well. 

We also discussed the general applications of the 
CPZM to high-Q resonance circuits. We showed that the 
response function of an arbitrary lumped or distributed 
resonator is or can be approximated by a ratio of poly- 
nomials. The ratio of polynomials can be approximated 
further to a linear factional form by finding its pole and 
zero closest to the bare resonance frequency of the res- 
onator. We described the general procedure for such ap- 
proximations. 

The methods presented here provide a systematic and 
rigorous treatment of response functions of high qual- 
ity factor resonators. They allow the reliable determi- 
nation of the quality factor and electric field in driven 
resonators. These methods are applicable to a variety of 
experimental investigations of the mechanisms of dielec- 
tric loss in the quantum regime. 



where N{lo, z) and D{lj, z) are polynomials of both lo 
and z. Provided (5w x r ^ 27r. we can take the first 



order expansion of z as z 



-h jreJ"o^(a; - w^). 



Thus we approximate the near-resonance response of a 
certain harmonic of the distributed resonator into a ratio 
of polynomials. Then, the same procedure to reduce the 
response function to the linear fractional form in Eq. (50) 
can be applied. As a result, the CPZM is applicable to 
distributed resonators as well. 
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Their coefBcicnts are defined as 



Appendix A 



This appendix contains the notation definitions and 
details of approximations of 521 and V^/VjjJ" in Sec. II 
and III. 

We first consider the response functions, 5*21 and V^/l^^ 
(Eq. (1) and (2)), of the circuit model shown in Fig. 1 b) 
and replace the circuit parameters in these expressions 
with the reduced parameters using Eq. (6). After per- 
forming parameter replacements and expression simplifi- 
cations in Mathematical^ , we obtain the following poly- 
nomial factional forms for 5*21 and v = V/V:^^: 



S21 



flQ + aiqx + a2{qxY 



bo + biqx + b2{qxY 4 
Co + ciqx 



b3{qxf + b4qx)^' 
+ C2iqx)^ 



do + diqx + d2{qxY -K dz{qxY + +di{qx) 



(Al) 
4- (A2) 



and 



ao = + 32(-ig + 7-a7)^ (A3) 

ai = + 32(l + 7-a7)^ (A4) 

02 = + 8(1 -h 7 - a-i)i (A5) 

Co = + 32^(7 - a7 + iai) (A6) 

ci= + 16e(-2(-l + a)7 + mC) (A7) 

C2 = - 8(-l + a)7e (AS) 



J 



60 = do = - 16i {iq-f + {i{2 + (-2 -t- a)a)-f + q{2 + a'^j- /3-f)) ^ + {ip[q + i-f) + a^(l - iq + 7)) C^) 

61 = di = + 8 (2i7 + AC + 2(2 + 2(-l + a)a - l3)-f£, - i {3a^l -f 7) - /3(2 + 37)) f) 

62 = ^2 = + 4 (3i7 + (2 + (2 - 2a + 6a^ - 5/3) 7) C - 3z {a^ - (3) {I + 7)^^) 

63 = ds = + 2i7 + 8 {a^ - /3) 7? - 2i {a^ - /3) (1 + 7)^^ 

64 = ^4 = + (a^ - /3) iC 



(A9) 
(AlO) 
(All) 
(A12) 
(A13) 



We then extract the coefficients of the zcroth and first 
order qx terms from the corresponding numerator and 
denominator of each response function, which allows us 
to rewrite the approximate response function in the linear 
factional form. We obtain: 



S'21 = r; X 



and -D = C X 



a4 


-JX 


6 + 


jx' 


a' 


+ jx 



JX 



(A14) 
(A15) 



where 



^ g- j(-l + a)7 
g(l -1-7 — 07) 
2j(-1 + q)7 + 2< 
2q{-l + a)j -jqaS,' 



(A16) 
(A17) 



and 



^ ^ ^, _ -2g7 + 2 (-(2 + (-2 + a)a)^ + jq (2 + a^j - (3^)) ^ + 2 {~f3{q + ^7) + a^q + j(l + 7))) ^ .^^g. 



q (-27 + 2j(2 + (2 + 2(-l + a)a - (3)j)( + (3^2(1 + 7) - /3(2 + 37)) e) 
-Ai{-1 + a)j^ - 2aS.'^ 
^ "-27 + 2z(2 + (2 + 2(-l + a)a - /3)7)e + (Sa^il + 7) - /3(2 + 37)) e ' 

I 

We rewrite the linear factional S21 and v by separating g from the other reduced parameters, which gives: 

A + j^+:ix 
021 = V 



(A19) 



C + f +j(i?+f)+jV 



(A20) 



and 



i = C 



— + j- + jx 



C+i+jiE + ^) 



jx 



(A21) 
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By Taylor expanding the coefBcients A, C, and E in 
terms of small parameters a and 7, we find that A fn 1, 
C w 1, and E ^ 0. We identify the resonance frequency 
as the minimum value of |S'2i|- By taking the first order 
derivative of 15*211 with A = 1, C = 1, and _E = 0, we 
find the extrema: 



-B^ + D^ + F^ + 2Dq ± ^(£>2 + {B- F)^)((5 - F)^ + (£> + 2g)^ 

2{B^F)q 



(A22) 



Only the extremum with "— " sign corresponds to the minimum position. As a result, the true resonance frequency 
of the circuit is 



Wo = CJg 1 + 



B^ +D^ + F^ + 2Dq - ^{D^ + {B~ FY){{B - Ff + {D + 2g)2 

2{B - F) 



(A23) 



r 



We take the resonance frequency as the resonance fre- 
quency without internal loss [q = 0) as wq = Wo(l — 
B/2). Therefore, the response functions S'21 and v can 
be rewrite in reduced forms in terms of the resonance 
frequency Wq and internal quality factor Qi = — = 



LUoRCo 



l-B/2 



'S'21, red — V 



and 



Wred = C 



1 + 2j 



Luoqi 



+ g.(l-i3/2) '^J q,{l-B/2) '^ -^J uoq. 



(A24) 



K I 3(./-B) 

9i(l--B/2) ^ 9.(1-5/2) 



2 - t^-CJo 



^ q,(l-S/2) +-^9,(1-3/2) ^"^ a;oq. 



(A25) 



In the end, we calculate the quantity of the reduced re- 
sponse function tired for uj = luq- We obtain: 



vo = C- 



9i(l--B/2) ^ qdi-B/2) 



1 + 



D 



J{F-B) 



(A26) 



gi(l-S/2) ' q,(l-B/2) 



The quantity of uq depends on the values of the parame- 
ters which appear in the expression of S'21. red and a linear 
scaling factor A. Wc define A as 



A-C 



K + jjJ-B) 
1 - B/2 



(A27) 



The full expression of A in terms of the reduced parame- 
ters is given below: 



A = 



47C(2ae-(-l + a)7(-2j + ae)) 



(-2 + (-1 + a)7) (-27 + 2j(2 + (2 + 2(-l + a)a - /3)7)C + (3a2(i 4. ^) _ ^(2 + 3^)) ^2) ■ 



(A28) 



r 



The calculations involving replacing parameters, simpli- 
fying expressions, extracting coefficients from polynomi- 
als, and expanding expressions into Taylor series in this 
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